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ABSTRACT 

We present a blind multi-detector multi-component spectral 
matching method for all sky observations of the cosmic mi- 
crowave background, working on the spherical harmonics ba- 
sis. The method allows to estimate on a set of observation 
maps the power spectra of various components present in CMB 
data, their contribution levels in each detector and the noise 
levels. The method accounts for the instrumental effect of 
beam convolution. We have implemented the method on all 
sky Planck simulations containing five components and white 
noise, including beam smoothing effects. 

1. INTRODUCTION 

The precise measurement of Cosmic Microwave Background 
(CMB) anisotropies is one of the main objectives of modern 
cosmology. The small temperature fluctuations of the CMB 
with respect to the pointing direction on the sky reflect the 
primordial density perturbations in the young Universe. The 
spatial power spectrum of those fluctuations depend on a set 
of important parameters, known as cosmological parameters. 
The accurate estimation of CMB power spectrum is thus of 
prime importance in cosmology. The Planck satellite (to be 
launched in 2007) will map whole sky emission with unprece- 
dented signal to noise ratio. 

The accuracy of CMB power spectrum measurements are lim- 
ited by the other astrophysical emissions present in the sub- 
millimeter range of the spectrum. Those emissions, called 
foregrounds (as they are emitted in front of the CMB), are of 
different origin. Some of them originate from within our own 
Galaxy, as the dust and synchrotron emissions, others are of 
extra-galactic origin as the Sunyaev-Zel'dovich effects. All 
depend on the wavelength. 



CMB from the separated CMB map. This approach is not 
fully satisfactory for two reasons. First, component separa- 
tion requires prior knowledge of the electromagnetic spectra 
of the components which are not all very well known. Second, 
it would be preferable to estimate the CMB power spectrum 
in one step by jointly analysing the observation maps. 
A new approach to process multi-detector multi-component 
data has been developed in papers [ 9 1 and iflOll . The method is 
based on likelihood maximization in the Whittle approxima- 
tion and takes additive noise into account (in the CMB obser- 
vations, a significant amount of noise is expected). The spec- 
tral diversity of the various components is used. The method 
has been implemented on the domain of small sky maps, in 
the flat sky approximation. 

In this paper, we have refined the multi-detector multi-compo- 
nent (MDMC) spectral matching method in several aspects to 
account for distinctive CMB observation features. First, as 
current observations cover as large a portion of the sky as 
possible, we have adapted the method for a spherical har- 
monic expansion of all sky maps in order to exploit all the 
information contained in the maps. Second, we account for 
the finite spatial resolution of the detectors, to take into ac- 
count the fact that the sky is not seen with the same resolution 
by all detectors. Finally, our method allows for the inclusion 
of some physical knowledge about the components. Fixing 
some parameters, the likelihood is maximized over the other 
parameters. This operation allows, for example, to break de- 
generacies due to components having similar spatial power 
spectra. 

We have implemented the method on full sky simulated Planck 
observations. We compare the estimated CMB power spectra 
in the blind approach and in the ideal case where one has a 
perfect knowledge of the component emission laws. 



It is thus important for a precise measurement of the CMB 
power spectrum to deal with the various foregrounds present 
with the CMB anisotropies. The availability of several detec- 
tors operating in several bands (10 for Planck ranging from 30 
to 850 GHz) allows to distinguish the various contributions. 
Component separation methods has been adressed by a num- 
ber of authors D]0l3HEl|6ll2IEI The standard approach 
consists in producing the cleanest maps of the various com- 
ponents, followed by estimation of the power spectrum of the 



2. MODEL OF SKY EMISSION 

The key assumption is that the sky emission at a given fre- 
quency is the linear superposition of astrophysical compo- 
nents. In addition, we assume that the emission laws of the 
various components do not depend on the position on the sky. 
The signal measured by a detector is the sky emission con- 
volved by a beam shape (depending on the detector), plus an 
additive noise. Assuming a symmetric beam shape, the obser- 



vation by a detector d is given by : 

x d(r) = Yl A dc / b d (\r-P\).s c (P)dP + n d (r) (1) 

c=l ^ 

r is the direction of observation on the sky, s c is the emission 
template for source c, rid represents the noise of the detec- 
tor cZ, bd (\r — represents the beam, and depends only on 
\f — r' | (or on r.r') for symmetric beams and A is the mixing 
matrix. Each element of the mixing matrix results from the 
integration of the emission law of one component over one 
detector frequency band. 

A natural basis for the application of the MDMC spectral 
matching method, exploiting the components spectral diver- 
sity is the spherical harmonics basis. We write the decompo- 
sition of the signal in this basis : 

oo I 

^(r) = J2 £ ^(?,m) W) (2) 

1=0 m=-l 

where Y l m (r) are the spherical harmonic functions. Coeffi- 
cients Xd(l, m) can be easily computed using the orthogonal- 
ity of spherical harmonic functions : 

x d (l,m)= I x d (r)Yr(?ydn (3) 

The convolution between a symmetric beam and the sig- 
nal in real space becomes a product in the spherical harmonic 
basis. Then, we obtain the observation coefficients combining 
equations [l] and |3 

x d (l, m) = b d (l) ^2 Adc Sc ^' m ) + m ) ^ 

c=l 

where bd (I) is the Legendre Polynomial expansion of 
bd(cosB), 6 being the angular distance from the center of 
the beam so that cos6 = r.r' . For a Gaussian beam, b(l) ~ 
exp(-of 1(1 + l)/2) and a b = 0beam/V81n2, where 6 b eam 
is the full width at half maximum of the beam. 
Let us consider the diagonal matrix B(l) such that the diag- 
onal element Bdd(l) — b d (l)- We define the new coefficients 
x'(l,m) = B(l)~ 1 x(l,m). It is useful to write these "decon- 
volved" observations coefficients in a matrix form : 

x'(l, m) = As(l, m) + B(l)~ l n(l, m) (5) 

The introduction of x' will be justified later. 
Spectral statistics 

We now need to compute the spectral statistics of the ob- 
servations x' (Z, m), given by R X >(1, m ) — {x\l,m)x'(l,m)^), 
where denotes transpose-conjugation. 

R x , (Z, m) = AC(l, m)A ] + BQ^NQ,, m)B(l)- ] (6) 

C(Z, m) and N(l, m) are respectively the component and the 
noise co variance matrices. Statistical independence between 
components implies that Ci m is a diagonal matrix. We also 



assume that the noise is white and independent across detec- 
tors, so that : TV (7, m) = diag(af J The main reason 
to work with variable x' is that we can average the covari- 
ance matrices R x >(l,m) over bins and preserve the simple 
structure of equation © for the signal part. Let us define the 
following particular averaging over bins: 

q l = lmin(q) m=~l 

Here q = 1, Q is the spectral bin index, the modes (Z, m) con- 
tributing to q are such that Z m i n (g) < Z < Z max (g). Since m 
can vary between — I and Z, the number of modes in each bin q 
is n q = Y^\Zif^(q) (2Z + 1). The reason for the choice of such 
averaging is that the spectral covariance matrices of isotropic 
components on the sky (we have strong reasons to think that 
cosmological components are isotropic) do not depend on the 
parameter m, so C(Z,m) = C(Z). The average covariance 
matrices are : 

R x >(q) = AC(q)Al+M{q) (8) 

where M(q) = l/n q Y^iZT^X) ( 2l + 1 ) B ( l )~ 2N is a dia g°" 
nal matrix. 

The band-averaged spectral covariance matrices are esti- 
mated by : 

^ Uax(g) I 

^to = — E x'(l,rn)x'(l,m)i (9) 

q l=l m in(q) m=-l 

which is real valued since x f (l,—m) = x f (l,my for real 
data. 

3. THE MDMC METHOD 

The aim of the MDMC spectral matching method is to obtain 
an estimate of different parameters of the model, of relevance 
in astrophysics and cosmology, without the need of prior in- 
formation. Those parameters are the mixing matrix A, the 
band averaged component power spectra C q and the noise co- 
variance N. They are collectively referred to as 
= {A,C(q),N}. 

The method is based on minimazing the mismatch between 
the empirical spectral covariance matrices of the observations 
(equation |9l> and their expected values which depends on the 
mixing parameters (equation 0. The mismatch is quantified 
by the average divergence measure between two matrices : 

Q 

where D(R\ , R2) is a mesure of the divergence (the Kullback 
divergence) between two positive matrices and n q is the num- 
ber of modes in each band q. Assuming that the spherical 
harmonics coefficients of the components are random real- 
izations of Gaussian field of variance R x > (q) and are uncorre- 
cted, the log-likelihood (up to irrelevant factor) take the same 



form as in equation [To] (in the frame of the Whittle approxi- 
mation). The divergence is given in this case by : 

D(R U R 2 ) =tr(R 1 R^ 1 ) - log detfT^i^ 1 ) - m d (11) 

The parameter estimate is given by 6 = argmin^ The 
connection with the log-likelihood guarantees good statistical 
properties of the estimates (at least asymptotically). 
Optimization algorithm: 

The optimization is made using an Expectation-Maximi- 
zation (EM) algorithm, and completed by quasi-Newton al- 
gorithm (BFGS) in order to accelerate the convergence. The 
formal algorithm described in 1 9 1 is used, with minor changes 
(the equations change slighly when we account for the beam 
smoothing effect). 

Degeneracies 

As seen from equation [8l one can exchange a factor be- 
tween the mixing matrix A and the component power spectra 
C(q) without changing the result of the equation. We fix this 
degeneracy in the EM step by fixing the norm of each column 
of A to unity; the power spectra are adjusted accordingly. In 
the quasi-Newton step, a penalty term is added to which 
penalizes the deviation of each column of A with a minimum 
penalization at unity. 

4. APPLICATION 

We now turn to the application of the MDMC spectral match- 
ing method on synthetic data. 

4.1. Simulated Planck data 

We use full sky Planck simulations provided by the Planck 
consortium. The maps are generated at the ten frequencies 
of Planck instruments (30, 44, 70, 100 GHZ for the low fre- 
quency instrument and 100, 143, 217, 353, 545, 857 GHz for 
the high frequency instrument). Five components and white 
noise at the nominal level of Planck instruments are mixed 
according to equation |4| (the beam sizes are by increasing or- 
der in frequency : 33, 23, 14, 10, 10.6, 7.4, 4.9, 4.5, 4.5, 4.5 
arcminutes). The components are the CMB, the kinetic and 
thermal SZ effect from galaxy clusters, and the galactic dust 
and synchrotron. They are obtained as follows: 
The CMB component is randomly generated using a power 
spectrum C\ predicted by CMBFast f\H with standard cos- 
mological parameters. The galactic components are gener- 
ated using observation templates at different frequencies, the 
galactic dust is modeled using the DIRBE-IRAS 100/im maps 
and the galactic synchrotron is a destriped version of the 408 
MHz Haslam survey with additional small scale structures 
(see (4J). The Thermal and the Kinetic templates are simu- 
lated with a gas dynamics code IT2l . Note that the kinetic SZ 
effect (always sub-dominant) and the CMB component have 
similar emission laws. The simulations are performed up to 
the resolution of 3.5 arcminutes. Figure shows the simu- 
lated Planck observation maps at all the frequencies. 



4.2. Estimated parameters 

We have applied our method on the above synthetic data. It 
is necessary to fix the number of components assumed to be 
present in the data. We choose to characterize 4 components. 
One of the components in the simulations, the kinetic SZ, is 
negligible at all frequencies. Moreover, it can not be separated 
from the CMB directly with our approach as the consequence 
of their proportional emission law (CMB and kinetic SZ form 
one component). 

Spherical harmonics are computed up to multipole /=3000 
(corresponding to 3.5 arcminutes), and we choose bins of 
width Al = 10. We investigate two different approaches : 

• First we adopt a quasi-blind approach. We estimate all the 
elements of the mixing matrix except three elements, corre- 
sponding to three of the four elements at 857GHz, which are 
fixed to zero (we expect the presence of only one component 
at this frequency). We will discuss the reason for this choice 
in subsection 14.41 All the other parameters are estimated, 
including component power spectra and the noise variances. 
The total number of parameters isl0x4-3 + 4x 300 + 10 = 
1247, compared to 300 x 10 x (10 + l)/2 = 16500 data ele- 
ments 

• In a second approach, we assume that we know the elec- 
tromagnetic spectra of the components perfectly, so we fix all 
the elements of the mixing matrix to their true values. We 
estimate the component spatial power spectra and the noise 
variances. 

4.3. Results 

In the blind case, all the parameters of the mixing matrix cor- 
responding to the CMB and dust components are estimated 
with a very good accuracy. Table [l] gives the ratio between 
the recovered and the true mixing elements of the CMB. The 
mixing matrix elements corresponding to the thermal SZ are 
estimated with a good precision. The galactic synchrotron 
emission law is very well constrained at lower frequencies. 
Therefore, we show that strong constraints can be put with 
our method on the component emission law, in particular for 
the galactic components at frequencies far from their maxi- 
mum emission. 

We now concentrate on the spatial power spectrum of the 
CMB component. Figure [l] shows the estimated CMB power 
spectrum in the blind approach. Figure shows the relative 
errors made on the CMB power spectrum estimation given by 
(\C(q) — C(q)\/C(q)) for the two approaches. The remark- 
able point is that the estimation errors in both approaches are 
equivalent. Therefore, it seems that the emission intensity of 
the components at the different observation frequencies can 
be estimated in Planck data without loss of accuracy in the 
linear CMB power spectrum estimation. 
In both cases, the method allows to estimate accurately the 
CMB power spectrum up to I ~ 2000. At smaller scales the 
dispersion begins to be significant. This result is not surpris- 
ing since the noise and the beam smoothing effect are impor- 
tant at these scales for all detectors. Also, the estimated power 



channel 30 44 70 lOO(LFI) lQQ(HFI) 143 217 353 545 

CMB 0.999984 1.000254 0.999780 1.QQQQ81 1 0.999993 0.999836 0.998972 0.990155 

Table 1. Ratio between the recovered and the true mixing elements of the CMB. The recovered parameters are rescaled to the 
detector at 100GHz of HFI (which is the best channel for the CMB). 
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Fig. 1. Estimated CMB power spectrum in the blind approach 



spectrum does not seem to be contaminated by the kinetic SZ. 



4.4. Discussion 

Semi-blind approach? 

As we have seen before, in the case where we estimate the 
mixing matrix as well as the power spectra of the components, 
we have fixed to zero the contribution of three components at 
857GHz, the only component we assume to be present at this 
frequency is the galactic dust. This hypothesis, being very 
realistic, appears to be necessary because we expect that the 
galactic dust and the galactic synchrotron have quasi propor- 
tional power spectra. Without any prior, the method can not 
separate these two components, but the previous simple pro- 
cedure of fixing some parameters allows to break the degen- 
eracies and to estimate accurately all the components. 

Partial coverage 

The galactic components are very inhomogeneous on the 
sky. In the region of the galactic plane, their intensities are 
several hundred times stronger than in the outer regions. Since 
we model the components as homogeneous, the method, when 
applied on those data, is sub-optimal. However, more ac- 
curate CMB power spectrum estimation may be obtained by 
making a galactic cut, involving a partially covered sky pro- 
cessing. The coefficients obtained from the spherical har- 
monic decomposition of partially covered sky are correlated. 
Our spectral matching method can be used on such data and 
a good choice would be to take bin size Al larger than the 
correlation length in the power spectrum measurement. The 
method, in particular, is applicable to data such that from the 
Archeops balloon-borne experiment in which the sky cover- 
age is about 30% of the whole sky. 




Fig. 2. Relative errors made on the CMB power spectrum 
estimation in the (quasi) blind approach (top) and in the semi- 
blind approach (bottom) 



5. CONCLUSION 



We have adapted our blind MDMC spectral matching method 
for the processing of all sky CMB maps. The observations are 
modeled as a noisy linear mixture of beam-convolved com- 
ponents. By maximizing the likelihood of the system, we 
estimate the power spectra of the components, their contri- 
bution levels at each frequency as well as the noise levels. 
The method has been applied on full sky Planck simulations 
containing five components. The spherical harmonics basis 
was used. 

The power spectrum of the CMB is accurately estimated up to 
/ c± 2000 in bins of size Al = 10. The comparison between 
the results obtained in blind and semi-blind approaches shows 
that the mixing elements can be estimated without loss of ac- 
curacy in the CMB power spectrum estimation. 




Fig. 



3. Simulated Planck observation maps at ten frequencies between 30 and 857 GHz 
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